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A  quantitative  study  is  made  of  the  relative  importance  of  the  purely  gravitationally  induced 
compression  (GIC)  due  to  fluid  particle  altitude  change  and  the  actual  “wave  compression"  which  can 
occur  at  a  fixed  altitude  in  a  gravity  wave.  The  results  for  density,  pressure,  and  temperature  variations 
show  the  following:  (1)  the  GIC  effects  predominate  (>95%)  for  vie  <20%,  where  v  is  the  horizontal 
phase  velocity  and  where  very  simple  formulas  can  be  obtained;  (2)  the  relative  importance  depends 
strongly  on  frequency  for  wave  periods  less  than  10  min  but  becomes  totally  independent  of  frequency 
for  periods  greater  than  20  min;  and  (3)  the  temperature  measurements  can  be  quickly  converted  to 
height  variations  wherever  the  GIC  effect  predominates;  in  general,  the  conversion  is  equivalent  to  the 
adiabatic  lapse  rate,  i.e.,  a  10°  temperature  variation  corresponds  to  a  height  change  of  1  km.  In 
addition,  the  total  kinetic  energy  density  can  be  simply  expressed  in  terms  of  height  variation  and, 
whenever  the  GIC  effects  predominate,  can  be  very  easily  obtained  from  temperature  measurements. 
An  interesting  by-product  has  been  that  for  waves  of  small  horizontal  phase  speed,  the  total  wave 
kinetic  energy  at  any  frequency  is  equal  to  the  kinetic  energy  of  the  natural  (Brunt)  oscillation  of  an 
air  parcel  with  the  same  vertical  displacement. 


> 


1 .  Introduction 

The  purpose  of  the  present  paper  is  to  determine  the 
degree  to  which  gravity  waves  (to  be  abbreviated  GW)  are 
subject  to  purely  gravitationally  induced  expansions  and 
compressions  (GIC)  from  altitude  changes  of  a  given  parcel 
of  air  and  the  degree  to  which  they  are  caused  by  purely 
“wave  compression.”  (By  purely  “wave  compression”  we 
mean  the  acoustic  component  of  the  gravity  waves.)  If  for 
certain  types  of  GW  the  GIC  effects  can  cause  most  (say 
90%)  of  the  GW  variations  so  that  we  can  neglect  “wave 
compression,”  then  the  resulting  simplification  [Dewan  et 
al.,  1988]  can  allow  us  to  use  very  simple  formulas  for  the 
GW  parameters  (density,  pressure,  and  temperature  fluctu¬ 
ations)  as  well  as  to  make  quick  estimates  of  other  GW 
parameters.  For  instance,  a  quick  conversion  from  temper¬ 
ature  fluctuations  to  height  fluctuations  can  be  readily  ef¬ 
fected,  and  the  latter  can  be  used  for  a  quick  estimate  of  the 
kinetic  energy  density.  This  is  also  very  useful  in  airglow 
measurements  where  interferometer  measurements  give 
temperature  directly  and  hence  airglow  height  changes  can 
be  estimated. 

We  will  begin  by  showing  that  the  density,  pressure,  and 
temperature  fluctuations  can  be  explicitly  expressed  as  a 
sum  of  two  terms:  one  describing  the  GIC  compression,  and 
the  other  the  “wave  motion.”  The  former  vanishes  in  the 
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absence  of  gravity  ( g  =  0),  while  the  latter  remains  in  altered  , 
form  and  provides  for  purely  longitudinal  acoustic  wave  k 
compression. 

For  g  ¥=  0,  the  relative  importance  of  these  two  effects 
depends  on  the  magnitude  of  the  horizontal  phase  velocity  of 
the  GW  as  compared  to  the  speed  of  sound.  For  small-scale 
GW,  where  by  definition  we  mean  the  phase  velocity  is  much 
less  than  the  speed  of  sound  and  the  wavelength  is  also 
small,  we  have  found  that  indeed  the  GIC  effects  predomi¬ 
nate.  Since  many  of  the  observed  mid-latitude  airglow  fluc¬ 
tuations  (including  OH)  involve  relatively  small  scale  GW 
[Witt,  1962;  Clairemidi  et  al.,  1985;  Taylor  and  Hapgood , 

1988]  with  low  horizontal  phase  velocity  v  <  70  m  s~' 

[Vincent  and  Reid ,  1983],  we  may,  for  most  cases,  assume 
that  the  density,  pressure,  and  temperature  variations  are 
GIC  and  a  quick  conversion  to  height  variations  may  be 
made.  Very  simple  formulas  for  the  kinetic  energy  density 
can  be  derived,  where  k}  >  0,  both  for  the  general  case  and 
even  more  so  in  the  GIC  approximation. 

We  should  point  out  that  there  is  considerable  interference 
between  the  two  effects  and  that  only  if  either  is  greater  in 
magnitude  than  the  other  by  a  large  factor  can  we  comfort¬ 
ably  neglect  the  smaller  terms.  To  test  the  limits  of  this 
criterion,  we  will  consider  three  specific  experimental  obser¬ 
vations:  the  first  being  the  airglow  observation  made  by 
Taylor  and  Hapgood  [1988]  (in  association  with  the  MAP- 
STAR  program)  which  easily  meets  our  criterion;  the  second 
being  the  wave  observations  by  Clairemidi  et  al.  [1985] 
which  also  meet  our  criterion;  and  the  third  being  the 
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observations  on  noctilucent  clouds  made  by  Witt  [1962] 
which  is  marginally  above  the  limits  of  our  criterion.  The 
results  show  that  for  all  three  sets  of  observations  the  simple 
GIC  model  is  surprisingly  accurate  and  that  one  may  safely 
assume  that  many  of  the  mid-altitude  density,  pressure,  and 
temperature  fluctuations  are  governed  by  GICs.  In  general, 
we  have  proved  that  in  the  region  v  «  C  where  the  GIC 
expansion  or  contraction  effects  dominate,  a  10°K  variation 
in  temperature  would  indicate  an  altitude  change  of  1  km  of 
a  parcel  of  air.  We  have  also  found  that  in  the  same  region 
the  kinetic  energy  density  ~  equivalent  to  that  of  the  free 
vertical  oscillation  of  an  air  parcel  with  the  same  vertically 
displaced  amplitude  and  is  independent  of  the  wave  fre¬ 
quency.  Thus  the  above  conversion  can  be  immediately  used 
to  calculate  the  kinetic  energy  density. 

2.  GIC  Expansion  and  Contraction 
When  a  fluid  element  is  displaced  by  a  vertical  distance  h 
in  a  hydrostatic  fluid,  its  density  is  compressed  by  an  amount 
Apa.  This  density  change  A pa  is  produced  by  a  change  in 
background  pressure,  A p  (where  A p  =  ( dp0ldz)h  =  c2Aps  and 
Apj  is  the  density  variation  due  purely  to  adiabatic  compres¬ 
sion),  minus  the  change  in  background  density, 
Apo  =  (dpo/dz)h.  One  can  show  that  in  a  uniform  isothermal 
background  atmosphere  with  scale  height  H, 

A pa  A ps  Ap0  (y  -  1\  h 

- = - = -  (1) 

P  o  p  o  P  o  \  y  J  H 

where  we  have  used  the  fact  that  c2p0  =  yp0  and 

1  <?Po  _  ,  _  1  dp o 

Po  dz  p0  dz 

3.  A  Comparison  of  the  GIC  and  Wave 
Compression 

In  a  GW  with  frequency  ui,h  =  wliw  (to  first  order)  where 
h  and  w  are  the  vertical  displacement  and  velocity  fields, 
respectively.  At  the  Brunt  frequency,  w  =  u>b,  the  atmo¬ 
sphere  would  just  oscillate  vertically  at  its  natural  frequency 
with  the  density  fluctuations  given  by  (1).  In  general,  by 
using  Hines’  [1960]  formula  for  Ap/p0,  h  =  wlito,  and 
horizontal  velocity  field  u  we  can  show  that  to  the  first  order 
(see  Appendix  1), 

«*> 

Po  \  y  J  H  \cj  \vj  Po  Po 

where  u  is  the  horizontal  particle  velocity,  v  (=vx)  is  the 
horizontal  phase  velocity,  and  A pc/p0  (={v/c)2  ( ulv ))  is  the 
density  fluctuation  due  to  wave  compression.  From  (2a)  it  is 
clear  that  Apa/po  — » 0  as  g  — *  0  since  H  =  c2yg,  and  the  terms 
Apc/p0  will  be  appropriately  changed  to  provide  for  the 
purely  longitudinal  wave  compression  at  all  frequencies  for 
sound  wave  propagation  in  the  absence  of  gravity.  As 
g  — *  0,  H  — »  «,  and  using  (A22)  we  obtain 

|Ap|  I  Apc|  |uv|  (M2  +  M2),a 


where  kx  is  the  horizontal  wave  vector  Equation  (26)  clearly 
shows  a  purely  longitudinal  wave  comp  asion  when  g  =  0. 
From  now  on  we  shall  only  be  concerned  with  g  *  0.  In  that 
case,  from  (2a)  the  relative  importance  of  Apa/po  and  Ap(/po 
depends  on  ( v/c )  and,  in  order  to  remain  well  below  the 
Lindzen  [1981]  limit  for  GW  saturation,  |«|  <  v  [see  Frills , 
1984].  Clearly,  for  small-scale  GW  where  v  <sr  C,  the  density 
variation  comes  mainly  from  GIC  compression,  Ap0/po-  On 
the  other  hand,  for  v  -*  C  and  close  to  the  limit  for  GW 
saturation  the  wave  compression  A pf/p0  is  large,  and  we 
need  to  include  both  effects. 

In  general,  since  both  Apa/po  and  A p,./p0  are  complex 
quantities  with  different  phases,  considerable  interference 
occurs  between  the  two  terms,  and  only  if  one  has  a 
magnitude  of  at  least  4-5  times  greater  than  the  other  can  we 
neglect  the  other  term. 

To  make  a  quantitative  comparison  between  A pa  and  Ap,  . 
we  plot  |Apa|/|Ap|  and  |Apr|/|Ap)  as  a  function  of  v/c.  It  fan  be 
shown  for  k 2  >  0  (see  Appendix  2)  that  >  ^ 

|Ap0|  [  (y-  D{1  -(v/r)2}  I1'2  ? 

- =  -  (3) 

|Ap|  [(y  -  1)  -  (<ula>b)2(vlc)2 

|Apcl  1  -  (io/wb)2 

|Ap|  (y  -  1)  -  (<i>la>h)2(vlc)2  \c/ 

Thus,  in  general,  |Apa|/|Ap|  begins  with  the  value  unity  (i.e., 
density  variations  are  due  entirely  to  GIC  compression)  at 
v/c  =  0,  while  |Apt.|/|Ap|  begins  at  zero  (i.e.,  no  wave 
compression).  As  can  be  seen  from  (3)  and  (4),  when  v/c 
increases,  the  behavior  of  |Apa|/|Ap|  and  |Ap,.)/|Ap|  depends 
entirely  on  w/a>b  for  any  value  of  v/c. 

For  long-period  gravity  waves  where  to  <K  aib,  we  may 
neglect  the  terms  involving  v)Uob  in  the  denominator  for  both 
|ApJ/|Ap|  and  |Apc)/|Ap|.  For  such  cases, 


j?-*0 


and  |ApJ/|Apj  drops  from  unity  as  v/c  increases  while  |Apr|/ 
|Ap|  increases  linearly  with  v/c. 

Figure  1  shows  this  behavior  for  gravity  waves  with 
periods  greater  than  20  min.  Actually,  there  is  little  variation 
with  period  from  20  min  on  up,  and  a  plot  at  r  =  120  min  falls 
exactly  on  top  of  the  corresponding  curves  at  20  min.  Both 
curves  are  plotted  from  (3)  and  (4),  but  there  would  have 
been  no  difference  if  we  had  used  (5)  and  (6).  At  r  =  10  min 
(Figure  2),  |Apr|/|Ap|  begins  to  deviate  from  a  straight  line  but 
only  slightly.  |Apa|/|Ap|  drops  somewhat  less  as  v/c  increases, 
showing  an  increase  in  importance  relative  to  wave  compres¬ 
sion  for  a  given  phase  velocity.  Furthermore,  the  crossover 
between  the  GIC  compression  and  the  wave  compression 
remains  fairly  constant  at  v/c  ~  55-60%.  In  general,  for  GW 
with  frequencies  in  this  range  (i.e.,  r>  10  min),  if  we  use  the 
approximate  criterion  that  |Apa(/|Apf|  has  to  be  greater  than 
a  factor  of  5.  then  v/c  has  to  be  less  than  15%  or  the 
horizontal  phase  velocity  has  to  be  less  than  45ms'1.  This 
covers  a  fairly  large  range  and  includes  most  of  the  observed 
mid-altitude  small-scale  GW  from  airglow  data. 

The  vertical  dashed  line  indicates  the  position  where  k 2  = 
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Fig.  1.  Plot  of  the  induced  density  variation  due  to  the  gravita¬ 
tional  (dash-dot-dot  curve)  and  the  acoustic  compressions  (solid 
curve)  as  a  function  of  v/c  for  wave  periods  of  >20  min. 


0.  It  is  easy  to  show  from  Hines’  dispersion  formula  that 
k*  =  0  corresponds  to 


To  the  left  of  the  vertical  line,  k*  >  0,  and  GW  may 
propagate  freely.  To  its  right,  k*  <  0,  and  the  GW  is 
evanescent  in  the  vertical  direction,  although  it  still  propa¬ 
gates  along  kx.  The  dotted  lines  represent  the  evanescent 
region.  At  v/c  =  1  we  have  the  Lamb  mode  for  which  there 
is  no  GIC  compression  since  there  is  no  vertical  motion  and 
variations  in  density  are  caused  entirely  by  wave  compres¬ 
sion,  an  exact  reversal  of  what  happens  at  v/c  =  0.  The  same 
reversal  occurs  in  all  the  figures. 

As  t-»  rb,  the  Brunt  period,  both  |Apa|/|Ap|  and  |Apr|/|Ap| 
increase  with  v/c  (see  Figures  3  and  4)  in  such  a  way  that  the 
GIC  compression  continues  to  predominate  over  a  much 
greater  range  of  horizontal  phase  velocity.  In  fact,  in  Figure 
4  where  r  -  5.2  min,  the  GIC  compression  continues  to 
dominate  all  the  way  to  the  dashed  line  where  v/c  ~  42%, 
corresponding  to  a  horizontal  phase  velocity  of  well  over  120 


Fig.  2.  Plot  of  the  induced  density  variation  due  to  the  gravita¬ 
tional  (dash-dot-dot  curve)  and  the  acoustic  compressions  (heavy 
solid  curve)  as  a  function  of  v/c  for  a  wave  period  of  10  min. 


m  s"1 :  a  strictly  medium-scale  GW  [Francis,  1975].  Thus  as 
the  GW  period  approaches  the  Brunt  period,  the  GIC 
compression  dominates  over  a  larger  range  of  v  and  hence  a 
larger  range  of  wavelengths. 

The  pressure  variation  can  be  easily  shown  to  be  the  same 
as  A pc.  Using  the  linearized  horizontal  momentum  conser¬ 
vation  equation,  it  can  be  shown  that  for  monochromatic 
waves, 

A  p  =  vp0u 
or 


so 

|AP|  jAPrl 
c2|Ap|  |Ap| 

From  Figures  1-4  we  see  that  there  is  little  pressure 
variation  produced  by  wave  compression  for  small-scale 
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Fig.  3.  Same  as  Figure  2  except  for  a  wave  period  of  6.4  min. 


GW.  Only  as  the  scale  size  of  the  GW  increases  does  it 
become  important.  Thus  for  small-scale  GW  (<45  m  s  ” 1 )  the 
pressure  variations  are  relatively  insignificant. 

The  temperature  variations  may  be  very  easily  obtained 
from  the  linearized  perfect  gas  law  given  by 


AT  Ap  Ap 
To  Po  Po 

From  (7),  (9),  and  (2a)  we  obtain 


AT  ( vu\ 


Thus  for  v  <k  c  (small-scale  GW), 


Apa 

Po 


(9) 


(10) 


AT 

To" 


Apa 

Po 


(y  -  I )h 

y  h 


(it) 


and  the  temperature  variation  is  opposite  in  phase  to  the  GIC 
density  compression.  We  will  define 


where 


AT a  Apa 
To  Po 


A  Tc  vu 

-  =  (7-1)- 
To  c 


Again,  as  in  (3)  and  (4)  it  is  possible  to  show  that  for  k2  >  0 
(see  Appendix  3), 


\ATa\ 

|A71 

jArc| 

|Ar| 


l-(v/c)2  ll/2 

1  -  (y  -  1  )((o/wb)2(v/c)2 

(y  -  1)(1  -  (Wan,)2)  m /v\ 
1  -  (y  -  l)(a>/<0 b)2(vlc)2  \cj 


(13) 

(14) 


Thus  once  again,  for  v/c  =  0,  |A7’a|/|A7’|  =  1  and  |A7'f.|/|7'|  = 
0,  and  we  see  that  in  the  low  phase  velocity  limit  (very  small 
scale  GW)  the  temperature  variation  is  caused  entirely  by 
the  adiabatic  expansion  and  compression  of  an  air  parcel 
oscillating  in  the  vertical  direction.  For  long-period  GW 
(w  <5C  <oh),  (13)  and  (14)  simplify  to 


AT_ATa  ATC 
To  To  T0 


\ATa\ 

|AT| 


21 1/2 


(12) 


(15) 
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Fig.  5.  Plot  of  the  induced  temperature  variation  due  to  the 
gravitational  (dash-dot-dot  curve)  and  the  acoustic  compressions 
(solid  curve)  as  a  function  of  vie  for  wave  periods  of  >30  min. 


Fig.  6.  Plot  of  the  induced  temperature  variation  due  to  the 
gravitational  (dash-dot-dot  curve)  and  the  acoustic  compressions 
(solid  curve)  as  a  function  of  vie  for  a  wave  period  of  10  min. 


|A71 


(16) 


Comparing  (15)  and  (16)  with  (5)  and  (6),  we  see  that  for 
long-period  GW,  |A7'a|/|A7"|  behaves  in  the  same  way  as 
|ApJ/|Ap|  in  (5)  while  |A7'f.|/|A7'|  has  a  far  more  gradual  slope 
than  |Apc|/|Ap|.  Thus  for  temperature  variations  the  GIC 
effect  predominates  over  the  “wave  effect”  through  a 
greater  range  of  phase  velocities. 

All  this  can  be  seen  from  Figure  5,  which  shows  the  behavior 
of  |A7„|/|AT|  and  |A7'C|/|A7'|  for  GW  with  periods  from  30  min 
on  up.  Again,  a  plot  of  both  curves  at  t  =  120  min  can  be 
exactly  superimposed  over  the  two  curves  for  r  =  30  min.  This 
is  similar  to  the  curve  of  the  density  variations  (Figure  1)  which 
has  a  lower  limit  at  t  =  20  min.  Again,  for  t  >  30  min  there  is 
no  significant  difference  between  (13)  and  (15)  or  (14)  and  (16). 

Figures  6,  7,  and  8  show  the  behavior  of  temperature 
variations  for  t  =  10, 6.4,  and  5.2  min,  respectively.  The  big 
difference  between  these  temperature  variations  and  the 
density  variations  given  by  Figures  2,  3,  and  4  is  that  for  the 
former  the  effect  of  GIC  compression  is  even  much  more 
important  than  the  wave  compression.  For  instance,  at  r  = 


6.4  min  the  wave  compression  effect  overtakes  the  GIC 
effect  at  v/c  =  0.72,  whereas  the  GIC  effect  remains  predom¬ 
inant  in  the  corresponding  temperature  variation  all  the  way 
to  the  vertical  dashed  line.  For  both,  however,  the  GIC 
effect  predominates  all  the  way  to  the  dashed  line  when  r  = 
5.2  min,  approaching  the  Brunt  period. 

The  dotted  curves  in  all  the  figures  on  the  right-hand  side 
of  the  dashed  line  correspond  to  evanescent  waves  with  k}  < 
0  and  thus  a  purely  imaginary  kz.  In  this  region  one  can  no 
longer  use  the  relatively  simple  expressions  given  by  (3),  (4), 
(13),  and  (14),  which  were  derived  for  kz  >  0  and  real  k,. 
Instead  one  has  to  work  with  the  general  polarization  rela¬ 
tions  given  by  Hines  [I960]. 

From  the  above  discussions,  it  is  clear  that  comparatively 
small  scale  GW  with  horizontal  phase  velocities  less  than 
40-50  m  s_1  may  be  considered  to  be  purely  under  the 
influence  of  GIC  expansion  and  compression  (>90%).  This 
will  allow  us  to  make  simple  conversions  from  the  measure¬ 
ment  of  temperature  fluctuations  to  height  fluctuations.  Thus 

|ATJ  |A7~n|  (y- 1)|*| 

T0  ~  T0  y  H  ( 
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Fig.  7.  Same  as  Figure  6  except  for  a  wave  period  of  6.4  min.  Fig.  8.  Same  as  Figure  6  except  for  a  wave  period  of  5.2  min. 


and  using  the  usual  values  for  y,  H,  and  r0,  we  obtain  the 
approximate  conversion  that  a  10°  variation  in  temperature 
corresponds  approximately  to  a  height  variation  of  1  km  at 
90-100  km  altitude. 

This  can  be  immediately  used  for  computing  the  time- 
averaged  kinetic  energy  density,  KE,  given  by 

KE  =  ip0(l«l2  +  M2)  (18) 


where  tv  is  the  vertical  particle  velocity.  Using  the  continuity 
equation  and  Hines'  dispersion  relations,  one  can  show  that 
(Appendix  4) 


KE  =  JPo 


1  -  (a>lwb)2(vlc)2 
1  -  (vie)2 


For  large-scale  GW  where  v  — >  c,  the  denominator  can  be 
much  less  than  the  numerator.  This  is  especially  true  for 
long-period  GW.  Thus  the  KE  can  then  be  much  greater  than 
^Ai«||/i|:.  For  small-scale  GW  (v  «  c),  on  the  other  hand, 
the  KE  is  very  accurately  given  by 

KE  =  \p0<ol\h\2  (20 a) 


It  is  interesting  to  note  that  (20a)  is  true  irrespective  of  the 
frequency  of  the  small-scale  GW.  Since  (20a)  is  also  the  KE 


for  the  natural  Brunt  oscillation  of  a  parcel  of  air,  we  may 
conclude  that  the  total  KE  of  small-scale  GW  at  any  fre¬ 
quency  is  equivalent  to  natural  vertical  atmospheric  oscilla¬ 
tions  with  the  same  vertical  displacement.  Here  the  GIC 
approximations  for  the  density,  pressure,  and  temperature 
variations  are  also  good.  Using  (17),  (20a)  acquires  the  form 


KE 

Po 


I  (  y  \  2  __  1  /  v  \  ^  2 

4  XT'  —  1/  T0  ~  4  y-y  -  1 J  T0 


(20  b) 


where  p0  is  the  background  pressure  energy.  Thus  a  mea¬ 
surement  of  the  rotational  temperature  variation  can  imme¬ 
diately  yield  the  total  kinetic  energy  density  at  any  frequency 
for  small-scale  GW.  In  fact,  (19)  shows  that  the  KE  depends 
strongly  on  frequency  only  for  large-scale  GW  where  v  -*  c 
(see  Figure  9). 

To  summarize:  (1)  for  small-scale  GW  (i>  <  45-50  m  s"1) 
the  KE  density  is  the  same  as  that  for  natural  atmospheric 
oscillations  and  is  independent  of  wave  frequency;  (2)  for 
any  GW  (arbitrary  v/c )  the  KE  again  becomes  equal  to 
natural  atmospheric  oscillations  when  w— *  <ob  as  expected; 
and  (3)  for  large-scale  GW  (v/c  — »  1)  the  KE  for  the  same 
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Fig.  9.  Plot  of  the  ratio  of  the  kinetic  energy  density  to  the 
background  pressure  as  a  function  of  vie  for  wave  periods  of  20,  10, 
6.4,  and  5.2  min. 


where  h  is  expressed  in  kilometers.  Using  the  GIC  approx¬ 
imation,  we  obtain 


l*r,| 

To 


(y  ~  1)  1*1 

y  H 


=  4.355  x  10_2|/i| 


(22) 


The  error  from  using  (22)  is  only  0.2%. 

The  second  set  of  data  were  taken  by  Clairemidi  et  al. 
[1985],  They  have  measured  waves  with  velocity  v  =  15.6  m 
s'1  and  a  period  of  48  min.  The  full  temperature  fluctuation 
without  approximation  is  given  by 


x 


10"2|/i| 


(23) 


Compared  with  the  GIC  approximation  given  by  (22),  we 
have  an  error  of  the  order  of  0. 1%. 

The  third  set  of  data  were  taken  by  Witt  [1962].  The  GW 
had  v  =  75  m  s_l  and  Ax  =  50  km,  and  the  rest  of  the 
parameters  are  the  same  as  the  first  data  set.  We  obtain 


|AT| 

- =  4.5  x  10“2|/i| 

To 


(24) 


The  GIC  approximation  gives  the  same  result  as  (22).  Thus 
the  error  for  this  example  is  3%.  The  somewhat  bigger  error 
is  understandable  because  at  75  m  s  " 1  the  horizontal  phase 
velocity  is  somewhat  above  the  45-50  ms-1  upper  limit  that 
we  mentioned  earlier  as  the  safe  limit  (i.e.,  vie  —  15-20%). 
Actually,  for  this  example  there  is  considerable  interference 
between  the  two  terms  on  the  right-hand  side  of  (10). 


5.  Conclusion 

The  present  paper  seems  to  show  that  for  small-scale  GW 
(say  vie  <  18%)  the  GIC  approximation  is  very  good  to 
above  95%  in  temperature  and  density  fluctuations.  The 
approximation  appears  to  be  applicable  to  many  of  the 
observed  mid-altitude  airglow  data.  Quick  estimates  of  and 
simple  formulas  for  GW  parameters  such  as  height  variations 
and  kinetic  energy  density  may  be  immediately  obtained 
from  measurement  of  temperature  variations. 


vertical  displacement  becomes  very  large  and,  unlike  the 
two  previous  cases,  can  also  become  very  strongly  fre¬ 
quency  dependent. 


Appendix  1 

We  can  derive  equation  (2a)  by  using  the  usual  expression 
for  Ap/p0  [Hines,  1960],  given  by 


4.  Applications 

To  specifically  apply  the  above  analyses,  we  apply  them  to 
three  specific  sets  of  observed  data:  two  sets  involve  GW 
parameters  well  within  the  above  mentioned  limits  for  the 
validity  of  the  purely  GIC  model;  the  other  involves  param¬ 
eters  above  the  borderline.  The  first  set  were  data  taken  by 
Taylor  and  Hapgood  [1988].  The  observed  data  include  a 
GW  with  a  horizontal  phase  velocity  of  u  =  21. 7  ms-1, a 
horizontal  wavelength  of  kx  -  26  km,  and  a  period  t  =  20 
min.  Using  g  =  9.8  m  s~2,  H  -  6.56  km,  y  =  1.4,  rb  =5.1 
min  (the  Brunt  period),  and  the  speed  of  sound  c  =  300  m 
sH,  we  find  that  if  we  employ  the  full  expression  given  by 
(10)  with  help  from  Hines'  [I960]  polarization  and  dispersion 
relations  we  obtain 


|AT| 


=  4.364  x 


10"2|A| 


Ap 

Po 


f  2, 

j  o)2kl  +  i 

(r  )8  x  2 c2  J] 

(Al) 


where /=  aez/2#,e**"-*  ‘ r). 

We  eliminate  a>  by  using  u>  =  kxv  and  substitute  in  for  the 
scale  height  H  =  c2/yg  to  obtain 


7.  *  {*'V’  +  *'V[( V)  5  '  i^3]}/  <A2) 

We  obtain,  by  adding  and  subtracting  (vlc)2(y  -  MyH)  in 
(A2), 

Hc2\y  2/ 


(21) 


/  (A3) 
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Introducing  the  vertical  displacement  held. 


H  itoH  x  \c/ 


(A4) 


|Apa| 

|Ap| 


(y  -  l)e 


so 


we  get 


,A5> 

Finally,  we  use 

u  =  kjc2v[kz  -  /(I  -  yl2)glc2]f 

to  obtain 


|Apa| 

|Ap| 


1  -  (<a/wb)2  +  (y  -  l)e  +  (y  -  2) 
(y  -  0(1  -  ( vie )2) 


1/2 


(y  -  l)  -  (to! to b)  (vie)2 


1/2 


(All) 


which  is  equation  (3). 

To  compute  |Apc|/|Ap|,  we  note  that 


Ap 

Po 


■/  (A  12) 


Ap  /y  -  l\  h  uv 
Po  \  y  /  H +  c2 


(A6) 


where,  again,  we  use  the  polarization  relation  for  u  [Hines, 
I960].  Substituting  for  tj, 


which  is  exactly  equation  (2a). 

Appendix  2 

From  now  on  we  consider  only  k2  >  0  (i.e.,  kz  real).  To 
obtain  equation  (3),  we  note  (using  (A4))  that 


APr/Po  =  vlc2[tokxc2(kz  +  it))]/ 


so 


|Ap< 


k;  +  v 


1/2 


(A13) 


(A  14) 


I/I 


'(271)(s)@")‘’Vl  (A71 

Using  (Al),  eliminating  kx  and  g,  and  substituting  in  t)  = 
-(1  -  y/2)glc2,  we  obtain 

zMWfc'y'W- 


|Ap|  [k2z  +  [((y  -  \)ly)(\IH)e  +  t,]2  j 
where  we  have  used  (A8)  and  (A  13).  Using  (A9),  we  obtain 


<*>l  ~  0)1 


2 - £  ~  T)2 


iApci  ftoi-oi2  r/ 

|AP|  “  1  c2 

+  5  +  ’)  1 * {("2  - 


+  [((y  -  l)/y)(l IH)]2e  +  2(y  -  \)V/yH)  ~  '}w 


«o2|/|  (A8)  so 


Using  the  dispersion  relation  [Hines,  1960]  and  letting  e  = 
c2lv 2  -  1,  we  have 


iApci  r  i  -  (oj/(ok)2  1 

|Ap|  {(y  -  1)  -  (tolto b)2(vlc) 2  j 


1/2 


(vie)  (A15) 


a>l  -  at2  o)l -w2  a)l -to2 

kz  2  2  2 

IT  C*  C 

where  t/2  =  (to2  -  <ol)lc2.  Thus 
|Ap 


e  -  |tj|2  (A9) 


y  ~  l\  L 
»l  “  v  y  )  H 

r-mi 


which  is  equation  (4). 


From  (10)  we  write 


Appendix  3 


2(y  -  1)tj 


-(^n 


A T  vu  (y  -  \)h  , 

—  =  (y  -  1)  — - — —  =  (y  -  l){i//c2[w*,c2(*;  +  117)] 

T0  c  y  H 

-  ik2g(vlc)2(c2lv2  -  1)}/  (A16) 

where  we  have  used  (A4)  and  (A13).  Thus 

A T/T0  =  (y  -  l)k2v2{kz  +  i[r)  -  glc2(c2lv2  -  1)]}/ 

and  since  A TaIT0  =  -  Apfl/po  (see  (12)),  we  obtain  by  using 
(A7) 


1/2 


(A10) 


|AT0| 


(y  -  \)(k2v2glc2)e 


Using  <ol  =  (y  -  l)g2/c2  and  dividing  it  from  the  numer¬ 
ator  and  the  denominator,  we  obtain 


|A7]  (y-x)ky[kl  +  (Tl-(glc2)£)2}''2 

Using  again  the  dispersion  relation  for  k2.  we  obtain 


(A17) 
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lAra!  ,  .i,  2l)g  g1  ]~l/2 


[1  -  ( u>2c2)lg 2  +  e],/2 


f  1  ~  (vie)2  V 
{l  -  (m2v2)/g2\ 


|Ara|  f _ 1. 

m  [l-(y- 


-  (v/c)2  ) 

l)(Ww„)2(t//c)2 


The  total  KE  then  becomes 


which  is  equation  (13). 

Equation  (14)  may  be  similarly  obtained.  The  denominator 
is  the  same  as  (13),  as  expected,  but  the  numerator  is  given 
by  |(y  -  l)uv/c2|,  which  is  identical  to  (4)  except  for  the 
additional  (y  -  1)  factor. 


kE-!p— a>p 

Obviously,  (A23)  is  equation  (19). 
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Appendix  4 

The  KE  density  is  given  by 

KE  =  Jp0(l"l2  +  M2)  (A19) 

From  the  work  by  Hines  [1960]  and  the  definition  of  rj,  we 
may  write 

u  =  <i)kxc2[kz  +  irj]/  (A20) 


w  =  -tok2c\\  -  (vlc)2]f 


|w|2  =  ^T^¥w: 


From  the  dispersion  relation 


/*A2  <Vb~<V2  /v\20)l~(O2 

W  ~~^2  \<7  "2 

hi2  / v\ 2  ~ 

k2  ~\c)  a,2 


Hence 


|m|2  +  Iwl2  =  \w\2<  1  + 


[(vfrw2)  -  1][1  -  (v/c) 
(1  -  (v/c)2)2 


'/c)2]j 
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